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We study equilibrium density and spin density profiles for a model of cold one-dimensional spin 
1/2 fermions interacting via inverse square interaction and exchange in an external harmonic trap. 
This model is the well-known spin-Calogero model (sCM) and its fully nonlinear collective field 
theory description is known. We extend the field theory description to the presence of an external 
harmonic trap and obtain analytic results for statics and dynamics of the system. For instance, 
we find how the equilibrium density profile changes upon tuning the interaction strength. The 
results we obtain for equilibrium configurations are very similar to the ones obtained recently by 
Ma and Yang^ for a model of fermions with short ranged interactions. Our main approximation is 
the neglect of the terms of higher order in spatial derivatives in equations of motion - gradientless 
approximation'^. Within this approximation the hydrodynamic equations of motion can be written 
as a set of decoupled forced Riemann-Hopf equations for the dressed Fermi momenta of the model. 
This enables us to write analytical solutions for the dynamics of spin and charge. We describe the 
time evolution of the charge density when an initial non-equilibrium profile is created by cooling the 
gas with an additional potential in place and then suddenly removing the potential. We present our 
results as a simple "single-particle" evolution in the phase-space reminiscing a similar description of 
the dynamics of non-interacting one-dimensional fermions. 
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I. INTRODUCTION 

The possibility of creating one and quasi-one-dimensional systems by confining cold atoms in cigar-shaped harmonic 
traps^Ti^ has raised interest in one-dimensional models of many body systems. These models have been intensively 
studied since 1970's. Standard perturbative methods developed in many-body theory are often not applicable to 
one-dimensional models because the low dimensionality effectively makes any interaction strong. On the other hand 
there are some non-perturbative methods which work specifically for particular (integrable) models in one dimension. 
The Bethe Ansatz approach, for example, was successfully used in constructing the complete thermodynamics of 
quantum integrable systems. However, this approach is not very suitable for studying the dynamics and correlation 
functions, due to the complexity of Bethe Ansatz solutions. To address dynamic questions the method of collective 
field theorySrii: was developed. It is essentially a hydrodynamic approach to many body systems. Many phenomena 
such as spin-charge dynamics^ , solitonsH, shock wavesi^-, spin density evolution*^-^ and sound wave propagation^ which 
have become of increasing experimental interest^"— '^"^i^^ can be studied via the collective field theory formalism. This 
formalism can be used both in combination with Bethe Ansatz for integrable models where it can even produce some 
exact results for dynamical problems or as a phenomenological approach where exact microscopic derivations are too 
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difficult or even impossible. The hydrodynamic approach allows to study truly nonlinear and nonperturbative dynamic 
behavior^ of many body systems. Linearized hydrodynamic equations are equivalent to the method of bosonization 
which was widely used for treating interacting systems in one dimensio n^^i^^ . 

In previous publications^'^- the collective field theory approach was applied to the well-known spin-Calogero 
model^^'^*' to study the coupled nonlinear dynamics of spin and charge. The model has a long range l/r^ inter- 
action and is not the easiest one to realize experimentally. On the other hand, the potential l/r^ should be considered 
as relatively short ranged in one dimension. The model has another advantage: it is integrable and the integrability 
is not destroyed by the presence of an external harmonic potential32, V{x) = ^-^^x'^. In contrast, for, say, the quantum 
integrable model of fermions with delta-interaction^^ the integrability is destroyed by an external harmonic potential. 

In this paper, we explore the effects of the harmonic trap on the collective behavior of the spin-Calogero model 
(sCM) at zero temperature. Its Hamiltonian is given by^^ j^°i^^i^^ 

Here and throughout the paper we take the mass of particles as unity and measure distances in units of an oscillator 
length I = yjh/muj and energy in units of huj. The operator Pij exchanges the positions of particles i and j^. 
The coupling parameter A is positive and N is the total number of particles. The last term in ([T]) is the harmonic 
trap potential. It prevents particles from escaping to infinity and corresponds to effective optical potentials used in 
experiments to keep particles. The above model ([1]) was shown to be integrable^. The fully nonlinear hydrodynamics 
for the above model ([T]) without an external trap has been investigated in Ref. S 

The paper is organized as follows. In Sec. [ll] we present the collective description of the microscopic model ([T]) 
in terms of collective fields and write equations of motion for these fields. Then we obtain static solutions: density 
and spin density profiles in Sec. IIIII We consider the dependence of these equilibrium profiles on coupling and find 
it to be very similar to the recent predictions of Ma and Yang^ for the model of fermions with contact interaction 
in a harmonic trap. In Sec. IIVI we show how hydrodynamic fields evolve when the system is perturbed from the 
equilibrium configuration. We model an initial non-equilibrium profile as the one obtained by cooling the gas with an 
additional potential, i.e., keeping what is commonly referred to in literature^ as a "knife" - in place (for examples of 
experiments involving "knife" see RefsHlandi). In this section (jlVp we solve the hydrodynamic equations of motion 
exactly. These equations written for dressed Fermi momenta are reduced to forced Riemann-Hopf equations. The 
solutions of these equations have a very simple form. It turns out that in the phase-space picture solutions are given 
just by a rotation by an angle t [tut in physical units) similar to a classical harmonic oscillator. We also study the 
density dynamics and see how an initial density perturbation evolves. The exact solution of the forced Riemann-Hopf 
equation is presented in Appendix |X] for reader's convenience. 



II. COLLECTIVE FIELD THEORY 



In this section we summarize the main results of the collective field theory for sCM^i^i^ following notations 
of Ref. 2. We also include an external harmonic potential which is done in a straightforward way. The microscopic 
Hamiltonian ([T]) is rewritten in terms of hydrodynamic fields: the density of particles with spin up/down pi.2 and their 
respective velocity fields vi,2- Although, it is possible to study the exact quantum hydrodynamics of sCM^!^, in this 
paper we neglect hydrodynamic terms with higher order of spatial gradients (gradientless approximation) and treat 
the equations of motion classically following Ref. '5. This description is referred to as a gradientless hydrodynamics 
and is applicable only for sufficiently smooth and slowly evolving field configurations, where terms with derivatives 
of fields can be neglected. A fully nonlinear gradientless hydrodynamics can describe nonlinear phenomena missed in 
conventional linear bosonization approach. In Ref. H , this theory was used to study the non- linear coupling between 
the spin and charge degrees of freedom;^ In Ref. [iSl the Emptiness Formation probability (a particular n-point 
correlation function) was calculated using instanton approach to the collective field theory of sCM. 

The collective field theory for sCM in gradientless approximation is remarkably simple and allows for separation of 
variables in terms of dressed Fermi momenta^ . Densities and velocities are expressed as linear combinations of dressed 
Fermi momenta. 

The aim of this paper is to extend this field theory by including an external potential, in particular, a harmonic 
potential. We find analytic solutions for both static profiles and dynamics of charge and spin densities of sCM 
in the presence of harmonic trap. The non-equilibrium initial configuration is realized by cooling the gas with an 
appropriate knife in place and then suddenly removing the knife as done in experiments (details are in Sec. IIIII and 
Sec. IIV[) . For simplicity, we focus here on a particular hydrodynamics sector of the model, i.e., we assume that the 
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following inequality is valid at any time everywhere in space 

\vi - W2I < tt{pi - P2)- (2) 

In addition, we assume throughout this paper that "1" labels the majority spin, i.e.. Mi > M2, where Mi 2 is the 
total number of particles with spin 1 and spin 2 respectively, i.e.. 



/ + 00 
Pi,2 dx. (3) 
-CXD 



We introduce particular linear combinations of these fields 

%,i,Li = vi±{X + l)7rpi ± A7rp2, (4) 
kB.2,L2 = (A + l)w2-Awi±(2A + l)7rp2, (5) 

referring to kji, l.i,2 as to dressed Fermi momenta (see Ref. 2 for details). For future convenience we invert (|4I5I) to 
get 

_ {km - kLi) _ KkR2 - kL2) 

' 2^(A + 1) 2^(A + 1)(2A + 1)' ^ ' 

khi + km 

V. = ^ , (7) 

km - kL2 

- M2ATT)' 

kL2 + km + Kkhi + km) 

= ^(aTT^ ■ 

In terms of dressed momenta in the sector ^ the hydrodynamic Hamiltonian in harmonic trap in gradientless 
approximation takes the form 

f +00 



^ " 127r(A + l) {^^1 ~ + 2ATT " } '^'^ 

dx. (10) 



1 f^°° 2 ( km - khi , km - ki 



2 I f^Rl - ft^Ll Km - KL2 

2j^^ lv2vr(A + l) 27r(A + l)(2A + l) 

The first term in (|10p was derived in Ref. 2 and the second term is due to the presence of the external harmonic trap. 
Notice, that the second term in ([TU)) can be written as \Pc dx where Pc = Pi+ P2 (|6l8p and represents the effect 
of the harmonic potential. 

The Poisson's brackets between hydrodynamic fields are given by 

{Pa{x),Va'{y)} ^ Saa'S'ix - y), (11) 

where tr, ct' = 1,2 are spin labels. 

Then ([TTjl along with (|4l5p give the following brackets for k's, (here a,/3 take values Rl, LI, R2, L2) 

{fca(x), fc/3(y)} = 2TrSaSal3S'{x - y) (12) 

with 

SB,l = -SLi = A + 1, (13) 

sm = -SL2 = (A + 1)(2A+1). (14) 

The collective Hamiltonian (fTU|) with Poisson's brackets ([T^ generate equations of motion for fields kt — {H, k} 
which turn out to be the forced Riemann-Hopf equations. Namely, for any k — km, fciL, fefl, fc2L we have 

kt + kkx = —X. (15) 

For the case of a more general external potential V{x) the right hand side of should be replaced by —d^V. One 
would expect to arrive at Riemann-Hopf equation^ modified by a force term due to the external potential V{x). It is 
remarkable, though, that the coupling constant A does not enter (jlSp at all. 
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III. STATIC SOLUTIONS 



In this section we consider static density and velocity profiles of the sCM in harmonic trap in gradientless approxi- 
mation. We give simple analytical expressions for these profiles for a system with an arbitrary spin polarization having 
a fixed number Mi of spin-up and M2 of spin-down particles. We describe static profiles in the form of phase-space 
diagrams. This description is very simple and gives a direct way to studying the dynamics of sCM (see Sec. IIVI) . 
Although these results are obtained for a particular ID model with long range interaction (sCM) the static profiles 
look very similar to the ones recently obtained by Ma and Yang for a fermionic model with contact interaction in 
harmonic trap>i. As mentioned earlier this similarity is, probably, due to the fact that the potential can be 
considered in one dimension as a relatively short ranged. 

To obtain equilibrium (static) density and velocity profiles we assume no time dependence in (|15|) and obtain, 

d^iP + x^) = (16) 

which is the equation of a "circle" in the phase space {x — k) plane 

k'^ + x'^ const. (17) 

The constant in p7|) depends on A, Mi and M2 and can be easily determined using (|6I8|) along with (jS]). We find 
that the dressed-momenta for spin-up (kRi^Li = ki) and spin-down (fc_R2,L2 = ^2) particles satisfy the equations of 
two circles respectively. 

kl+x^ = 2(A + l)Afi + 2AAf2, (18) 
kl+x^ = 2(2A + 1)M2. (19) 

The equation (jlSp defines a double- valued function ki(x). The positive value is identified with fcfli(a;) while the 
negative one gives kLi{x). The values of kii2,L2{x) are obtained in a similar manner from the second equation (|19p. 

In order to write the expressions for charge-density (pc = pi + P2) and spin-density (ps = pi — P2) it is convenient 
to use re-scaled coordinate ry and re-scaled charge and spin densities pc,s given by 

(20) 



V(2A + l)iV' 

/(2A + 1) , , 

Pc,s = ^ ' Pes- (21) 

Here the total spin M, the total number of particles (charge) N and the magnetization u are defined respectively as 

M = Ml - M2, (22) 
N = Ml + M2, (23) 



M 

77' 



(24) 



From (|22r24p we have 



Mi,2 = ^(l±i^). (25) 
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With these definitions we have from (|18ll9p and ([6 



(2A + 1) I V : 1 , , , 



Ps 



(2A + 1) / r V 



7r(A + l) VV 2A + 1 

One could also notice as a cross-check that 

+00 2^ /'+00 



+00 



1 + — — -7?2_yT^7^ . (27) 



Pcrf?7 = ^ / Pcdx^ 1, (28) 
Psdrj ^ — I psdx ^ V. (29) 
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Figure 1: Phase-space picture for sCM with A = 2 in equihbrium with an overall magnetization v — 0.8. The radii of circles 
are given by y^N{2X + 1)(1 - ty) and y^N{2X + 1 + ly) for inner and outer circles respectively. Particles fill the phase space 
uniformly with the density of two particles (up and down) per the area 27r(2A + 1) in the inner circle and of one particle (up) 
per the area 27r(A + 1) in the annulus area between inner and outer circles. 

The analytical expressions for momenta ()18|19p and the analytical expressions for charge density ([26]) and spin 
density ([?7)) are the main results of this section. We see that in terms of phase space (Fig. [IJ the system can be 
described by two circles of different radii. The radii depend on the coupling strength A and the numbers of particles 
Ml and M2 as can be seen from (|18I19|) . This is consistent with an exclusion statistics picture'^ so that the particle 
with spin up occupies the area 27r(A + 1) in the phase space while in the domain where both spin up and spin down 
particles are present the area is 27r(2A + 1) per 2 particles. In the limit A = (free fermions) the inner circle in Fig. [T] 
is filled with double density compared to the annulus region between the inner and outer circles. This is reflected as 
a bump feature for the charge density (see Fig. [2]). In the limit of a very strong repulsion A — ?> +00 the particles are 
mutually exclusive with the phase-space area approximately 27rA per particle and the charge density does not show 
any bump (see Fig. [5]). For an arbitrary A the bump feature interpolates between these two limits as shown in Fig. [2] 
with an explicit formula for the static charge density profile given in Eq. (|26p . We notice here that the size of the 
cloud Lcigar — 2^jN{2X -I- 1 -I- ^) given by the support of (^5)) grows with A (ID gas expands when A increases). We 
have taken this main dependence into account by plotting density as a function of re-scaled coordinate 77. 

It is interesting to note that although the bump in charge density profile disappears in the limit of a very strong 
repulsion, the analogous feature in the spin-density profile is present for any A. It is practically intact when plotted 
in terms of re-scaled variables (see Fig. [3]) and is given explicitly by ([77)1 . 

One should expect qualitatively similar profiles for the fermions with short range repulsion in a harmonic trap. 
Indeed, the physical origin of the bump feature is transparent for noninteracting fermions. It comes just from a 
superposition of densities of clouds of different size. We expect that the repulsive interaction smears the bump feature 
making sure that the particles of different species avoid each other similarly to the Pauli exclusion of particles of the 
same species. Indeed, recent calculations of charge density profiles for fermions with contact interactions by Ma and 
Yang^ give results very similar (qualitatively) to the ones shown in Fig. [21 It would be interesting to check whether 
the dip in the spin density profile does not disappear with an increase of the repulsion strength in the model of Ref. 

a 

The following comment is in order. As the origin of the bump feature in the equilibrium charge density profile can be 
traced to the model of noninteracting fermions, this feature is generic and should be observed also in three-dimensional 
harmonic traps which are widely used in cold atom experiments. Indeed, the superposition of two ellipsoid-like clouds 
of spin-up and spin-down particles will give the bump in the overall number density of atoms. The repulsive interaction 
will smear the feature while the attractive one will amplify it. 

To manipulate cold atom systems experimentally additional external potentials are often used. The gas is cooled in 
the presence of a harmonic trap and an additional external potential. The latter is usually created by optical means 
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Figure 2: Equilibrium charge density profile for various values of coupling constant A for fixed magnetization v — 0.8 is shown 
in re-scaled variables (|20I21[I . A = corresponds to noninteracting fermions. Upon increasing the interaction strength A the 
equilibrium profile eventually loses its bump feature. This prediction for sCM is very similar to recent predictions of Ma and 
Yang for fermions with contact interaction^. 




Figure 3: Equilibrium spin density profile for various values of coupling constant A for the fixed magnetization v — 0.8. In 
re-scaled variables the dip in the spin density profile depends weakly on the interaction strength A. Compare with the charge 
density profile of Fig. (2] 



and referred to as "knife". It is possible to create an external potential different for different particle species (spin up 
and spin down here) . In the following we concentrate on the potential acting only on "charge" degrees of freedom. 
We choose a knife potential of the form Vknife — —Vqc^^ . Then the static equation Pd|) is modified as 

9J Y + Y+14m/ej =0. (30) 

Due to the inclusion of the knife (|30)) the circles in phase-space acquire peaks (attractive knife Vq > 0) or dips 
(repulsive knife Vq < 0). There is also a peak or dip in the corresponding charge density profile. The knife shape, 
the corresponding phase-space and the charge density profiles are shown in Fig. 21 After using the knife to create 
an initial density profile with a peak (dip) one can remove the knife potential and study the evolution of density (or 
velocity) profiles as a function of time. The initial configuration is a non-equilibrium configuration. Its dynamics is 
governed by ()15|4|5p . We study the corresponding evolution in the next section Sec. IIVI 



IV. DYNAMICS 



Once the knife potential (Fig. HJ is suddenly removed we expect the system to evolve. In this section we study 
this dynamical behavior. Usually, the hydrodynamic equations are coupled partial differential equations. For sCM 
in gradientless approximation it is possible to separate variables using dressed Fermi momenta instead of charge and 
spin densities and velocities^ . The dressed momenta are related to hydrodynamic densities and velocities by (|4l5p and 
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Vknife=-Vo^'-^V«^ 






Figure 4: (left to right) a) Attractive knife potential in the presence of which the fermi gas is cooled, b) The distortion of phase 
space circles due to the knife, c) The corresponding charge density profile. The values of v and A are v = 0.5 and A = 2. 



satisfy the simple forced Riemann-Hopf equation p^ . We reproduce it in this section for reader's convenience. 

kt + kk^ = -X, (31) 

where k = km, kiL, k2R, k2L- 

It is remarkable that the force term (right hand side) is the same for all dressed momenta. This is just another 
manifestation of the noninteracting nature of Calogero-type models (free particles with exclusion statistics). Of course, 
physical quantities of interest such as densities and velocities are some A-dependent linear combinations of k' s (|6I7I8I9P 
and hence do depend on the coupling strength A. 

As it is explained in Appendix the forced Riemann-Hopf equation (|3ip can be easily solved in parametric form. 
Given an initial profile k^^x) = k(x, t — 0) (specifying the curve in a phase space picture) we can write a profile k{x, t) 
at time t in a parametric form 



x{s;t) = i?(s) sin + Q!(s)] , 
k{s;t) = R{s) cos [t + a{s)] . 

Here the parameter is s and the functions i?(s) and a{s) are determined by an initial profile ko(x) as 



a{s) — tan ^ 



ko{s)J ' 



(32) 
(33) 



(34) 
(35) 



consistent with p2|33|) at i = 0. 

We immediately notice that the evolution p2l33p is just a rotation of the curve k{x) in the x — k phase space with 
constant angular velocity 1 (cj in dimensionfull variables). 

In this section we take the initial profile fco(a;) as the one obtained as an equilibrium profile in phase space in 
the presence of the knife potential (Fig. HJd). When the knife is removed suddenly this profile serves as an initial 
non-equilibrium profile. The time evolution in the phase-space is just a rotation of an initial profile by an angle 
6 = t as shown in Fig. [S] From the phase space picture at time t obtained by this rotation of an initial profile it is 
straightforward to compute the charge density evolution using (and similarly for spin and velocities evolutions) . 

The top row of Fig. O shows the phase space rotation at various times. The bottom row is the extracted charge 
density. Since we are dealiirg with a gradientless theory we can study profile evolution only at times for which the 
field profiles are smooth. This gradientless approximation is commonly employed in studying nonlinear equations^ 
and allows to study the evolution for a finite time when the nonlinear terms dominate the terms with higher order 
in spatial gradients (dispersive terms). Of course, this is possible only if an initial profile is sufficiently smooth. For 
longer times, the solution inevitably evolves towards configurations with large field gradients (such as shock wavesi^iii) 
and the gradientless approximation becomes inapplicable. Choosing a sufficiently broad profile of the knife potential 
we make sure that duriirg the initial stage of the evolution, corrections due to gradient terms in equations of motioir 
are small. We emphasize that this evolution obtained in gradientless approximation already shows some interesting 
features. 
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Figure 5: Top Row: (left to right) Evolution of phase space for time t=0,0.08 and 0.23 respectively. We see that this is merely a 
rotation by angle t. Bottow Row: Corresponding charge density evolution for times t=0, 0.08, 0.23. The additional peak created 
by the attractive knife flattens and eventually splits into two peaks. The values of v and A are v — 0.5 and A = 2. 

We see that the central peak (created by cooling with attractive knife) in the charge density profile slowly fiattens 
and eventually splits into two peaks (see the bottom row of Fig. [5]) . Due to nonlinear effects these two split peaks 
start to steepen and we expect that at that point, gradient corrections will play a role and the profile would develop 
dispersive shock waves. 

Although the solution of ([5T|) is not well defined beyond some time (gradient catastrophe time) the parametric 
solution (j32l33[) can be formally extended beyond that time and produces multiply-valued solutions. These multiply- 
valued solutions should not be used as the equation (PT|) has corrections with higher power of gradients which will 
significantly change the solution beyond the gradient catastrophe time. It is interesting however that in time 1 = 2?: 
the solution should reproduce an initial profile and this is an exact feature of the spin-Calogero model. Of course, any 
corrections to the sCM model destroying integrability will lead to equilibration of the system and the time-periodicity 
will be lost. Probing the charge dynamics in fermi gases at large times is indeed possible experimentally^ and it would 
be interesting to see this equilibration experimentally. 

While our analytical expressions p2j|33p are correct for free fermions and sCM (long ranged interactions) we expect 
qualitatively similar dynamics at small times for systems with short range interaction. The corresponding dynamics 
studies for fermions with contact interaction (short ranged) could be an interesting extension to the recent work of 
Ma and Yang=. To the best of our knowledge there are no such dynamic studies for fermions with contact interaction 
in external harmonic trap. We would also like to remark that although our predictions for dynamics are for pure ID 
systems we do expect a qualitatively similar behavior for the more familiar quasi ID (cigars) experimental setups^"—. 
We leave a quantitative analysis of quasi-one dimensional systems for future work. 



V. 



CONCLUSIONS 



In this paper we addressed statics and dynamics of a model of one-dimensional spin 1/2 fermions interacting 
through a long range inverse square interaction in an external harmonic trap ([T]). While one-dimensional systems 
with long range interactions are yet to be realized there has been a great progress iir experimental studies^"— of 
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quasi-one-dimensional fermionic models with contact-like interactions in an external harmonic trap. 

Inclusion of an external harmonic trap potential is known to break integrability of most models (for example models 
of bosons or fermions with delta interaction). On the contrary, the spin-Calogero model (sCM) ([T]) remains integrable 
even in the presence of an external harmonic potentiaP^. Similarly, the collective field theory of sCM with harmonic 
trap retains a rather simple structure with dynamics analogous to the one of non-interacting fermions in harmonic 
potential. In this paper we used the spin-Calogero model as a toy model of more realistic systems of cold Fermi atoms 
in quasi-one-dimensional harmonic traps. 

Using the collective field theory reviewed in Sec.|lT]we studied static density profiles in Sec. Illll as well as dynamics 
of the model (see Sec. IIV|) . 

The obtained static solutions are found to be qualitatively similar to the static solutions for fermions with contact 
interaction in harmonic trapi. This similarity with short range interactions, probably, can be attributed to low- 
dimensionality where can be considered as relatively short ranged. The obtained density profiles are given in 
(|26l27p and are shown in Figures [U [2] and [S] The solution of a similar problem for fermions with contact interactions 
requires numerical solution of Bethe Ansatz equations in combination with the Thomas-Fermi approximation^. An 
interesting feature of an equilibrium solution is the "bump" in the charge density profile present for a system with 
non-zero polarization. The reason for this feature is relatively straightforward as it is already present in the system of 
non-interacting fermions. It is interesting, however, that the smearing of this bump with an increase of the interaction 
is qualitatively very similar for both fermions with contact interaction^ and for sCM model (Fig. [5]). We also notice 
that while the bump feature in charge density (Fig. ^ eventually disappears at strong coupling the spin density 
profile remains robust (Fig. [3]) . It would be interesting to have a calculation of the spin density profile for the contact 
interaction model considered in Ref. [Ij. 

To study the dynamics of the cold gas in harmonic trap we create an initial non-equilibrium density profile by 
"cooling" the gas in an additional attractive potential ("knife"). We choose the form of the attractive knife potential 
Vknife = — Voe~^ similar to the one used in experiments-. When the knife is suddenly removed the density profile 
shown in Fig. |3] serves as an initial non-equilibrium profile which is expected to evolve in time. We show in Sec. IIVI 
that the central peak in the charge density profile (created by knife) slowly fiattens/broadens and eventually splits 
into two peaks (see the bottom row of Fig. [5]). Due to nonlinear effects these two split peaks start to steepen and 
we expect that at that point, gradient corrections will play a role and the profile would develop dispersive shock 
wavesi^iii. 

In conclusion, we studied equilibrium configurations as well as dynamics of spin-Calogero model in harmonic trap. 
We argued that the model can serve as a toy model for cold Fermi atoms in one dimensional traps due to the relatively 
short range nature of an inverse square potential in one dimension. The integrability of spin-Calogero model is not 
destroyed by the presence of harmonic potential and simple analytic solutions of hydrodynamic equations of motion 
for this model in gradientless approximation are readily available. 
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Appendix A: Exact solution for Riemann-Hopf equation in arbitrary potential 

In this appendix we find an exact solution k{x,t) of the forced Riemann-Hopf equation, 

kt + e'{k)K = -V'{x). (Al) 

Here e(fc) is some function referred to as "dispersion" and the prime means the derivative with respect to the indepen- 
dent variable (e.g., e'(fc) = de/dk). The problem is to solve (|A1[) with an initial condition k = fco(x) at i = 0. This 
problem can be easily solved using the general method of characteristics^^. In this appendix we solve the Euler-type 
Eq. (jAip rewriting it in the Lagrange formulation. Namely, we notice that the left hand side of (|A1[) can be interpreted 
as the time derivative in the reference frame of the moving particle. In other words (|Aip is equivalent to 



X = e'(fc), 
k = -V'{x). 



(A2) 
(A3) 
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Here dot means time derivative and k{x{t),t) = dtk + idxk is nothing else but the left hand side of (jAl[) . 
The system of equations (jA2IA3[) has a first integral of motion (energy) 

e{k) + V{x) = E, (A4) 



where E is a time-independent constant determined by an initial condition. We invert (jA4ll and use it and ()A2p to 
find a solution of (|A2IA3I) and, therefore, of (Ml 

k{x) ^ e-^{E-V{x)), (A5) 

E = e(fco(s)) + F(s), (A7) 

where the last equation gives the value of the energy E in terms of initial data parametrically given as a; = s and 
k = ko{s). 

The system (|A5IIA7[) defines the solution k{x,t) of (jAll) for a given initial profile kQ{x). 
Let us now assume that there is no external potential V{x) — and the dispersion is quadratic e{k) = 
Excluding E and s from (jA5IIA7P we have 

k^ka{x~kt), (A8) 

which is being solved with respect to k gives a well-known exact solution k{x,t) of (|A1[) . The solution 
extensively used in Ref. ,2: for studies of the dynamics of sCM in the absence of an external potential. 

In the presence of an external harmonic potential V{x) — uj'^x'^ /2 (and e(fc) = fc^/2) we have from (|A5IIA7|) 



k{x) = y/2E-uj^x'^, (A9) 

t = r , (Alo) 

Js ^2E - w2y2 

^ ^\ [fco(s)'+^V] . (All) 

Calculating the integral (jAlOp and excluding E from the system (jAQIIAllj) we obtain after straightforward manipula- 
tions 

ujx = i?(s) sin [wt + a(s)] , (A12) 

k = i?(s) cos [cjf + q;(s)] , (A13) 

where we introduced 



i?(s) = V^252 + fco(s)2. (A15) 

The equations (|A12IA13p together with definitions (IA14IA15|) give a parametric (s is the running parameter) 
solution of (|A1[) for the case of quadratic dispersion and harmonic external potential. Putting w 1 in (jA12IIA15|) 
we reproduce p2j|35p used in the main body of the paper. 
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